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The last few years have seen significant experimental progress in characterizing the 
copper-based hole-doped high temperature superconductors in the regime of low hole 
density, p. Quantum oscillations, NMR, X-ray, and STM experiments have shed much 
light on the nature of the ordering at low temperatures. We review evidence that the 
order parameter in the non-Lanthanum-based cuprates is a d-form factor density-wave. 
This novel order acts as an unexpected window into the electronic structure of the pseu¬ 
dogap phase at higher temperatures in zero field: we argue in favor of a ‘fractionalized 
Fermi liquid’ (FL*) with 4 pockets of spin S = 1/2, charge +e fermions enclosing an 
area specified by p. 
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1. Introduction 

Ever since their discovery 1 in 1986, the cuprate high-temperature superconductors 
(SC) have been among the most well-studied correlated electron materials, both 
theoretically and experimentally. These materials show a whole zoo of different 
phases as a function of hole doping. The parent antiferromagnetic Mott insulator, 
the d-wave superconductor and the Fermi liquid state at large doping are reason¬ 
ably well understood at this point. On the other hand, the phases that continue 
to elude an explanation, and occupy a large portion of the phase diagram (Fig. 1), 
are the ‘pseudogap’ and the ‘strange-metal’. It is essential to develop a microscopic 
understanding of these phases in order to resolve the mystery of cuprate super¬ 
conductivity, as they are the ‘normal’ states out of which superconductivity (and 
various other ordering tendencies) arise at low temperatures. 

Although we don’t have a complete understanding of these phases, a lot is now 
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Fig. 1. A temperature (T) vs. hole-doping (p) phase diagram for the non-La-based cuprates. The 
parent state is an antiferromagnetic (AF) Mott insulator which upon being doped leads to a d- wave 
superconductor (SC) with transition temperature, T c . The extremely overdoped phase is a familiar 
Fermi-liquid (FL). The two mysterious phases (denoted ‘?’) are the pseudogap (PG), which onsets 
at a high temperature (T*), and the strange-metal (SM), occupying a wide fan-shaped region. 
Multiple recent experiments have investigated the nature of the lower (Tpw < T*) temperature 
density-wave (DW) phase in the underdoped cuprates, which will be a central topic of this review. 


known phenomenologically about the pseudogap regime of the underdoped cuprates. 
The first signatures of the ‘high’-temperature pseudogap (T < T*) came from mea¬ 
surements of the Knight-shift 2 , followed by other spectroscopic measurements 3 , 
which showed evidence for a suppression in the density of states at the Fermi-level 
along specific regions of the Brillouin zone. However, there is no clear evidence for 
any broken symmetry in this regime. On the other hand, in recent years many re¬ 
markable experiments have revealed the nature of the ‘low’-temperature pseudogap 
(T < T** < T*), where a number of broken symmetries have, in fact, been observed. 
A major breakthrough was the initial observation of quantum oscillations at high 
magnetic-fields and very low temperatures, followed by the subsequent discovery of 
short-ranged charge-order even at zero field in the underdoped regime. This was 
marked by initial optimism that, perhaps, the “hidden” broken-symmetry phase 
responsible for giving rise to the pseudogap has finally been identified. However, 
a more systematic study has now made it quite clear that the somewhat fragile 
charge density wave itself can’t be responsible for causing the pseudogap; but it 
possibly arises at low temperatures out of the normal pseudogap phase with no 
broken symmetry (there may be broken discrete symmetries, such as time-reversal 
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or nematicity 4 , at higher temperatures, but these have little direct effect on the 
electronic structure or spin excitations). One might then ask, what have we really 
learnt post the discovery of charge order? Could it perhaps serve as a “window” 
into the nature of the pseudogap state itself? 

In this review, we shall take the point of view that, indeed, the micro-structure of 
the charge-order, namely its ‘form-factor’ and (doping-dependent) wavevector, shall 
provide us important clues about the nature of the higher temperature pseudogap 
phase. We will investigate low-energy particle-hole instabilities of possible normal 
states, and this will help us identify the necessary features required to give rise to 
a state that is qualitatively similar to the one observed in the experiments. This 
will also aid us in distinguishing between some of the different approaches that have 
been used to describe the pseudogap phase. 

Broadly speaking, theoretical approaches for characterizing the pseudogap can 
be classified into two major categories. In the first category, the build up of antifer¬ 
romagnetic correlations and the opening of a spin-gap 2 at T* signal the onset of a 
quantum spin liquid 5 9 . Moreover, as a function of decreasing temperature, there 
could be multiple crossovers within the metallic spin-liquid and there could also be 
instabilities to other symmetry-broken phases a at lower temperatures. However, in 
order for this to be a useful and precise characterization, there should be remnants 
of the ‘topological order’ of the spin liquid at high temperatures 13 . One possibility 
is the presence of closed Fermi pockets which violate the Luttinger theorem con¬ 
straining the total area enclosed by the Fermi surface 10 , and this may be related 
to photoemission spectra which have intensity only on open arcs in the Brillouin 
zone. In this approach, the pseudogap metal found at high temperatures is a novel 
quantum state which, with moderate changes, could be stable at low temperatures 
for suitable model Hamiltonians. 

In the second category 11-14 , the antiferromagnetic correlations are precursors to 
the appearance of antiferromagnetism, superconductivity, charge density-wave, and 
possibly other conventional orders at low temperatures. In the pseudogap regime, 
we then have primarily thermal and classical, rather than quantum, fluctuations of 
these orders. This raises an immediate question: why don’t other unconventional 
superconductors, many of which have a robust antiferromagnetic phase 15 , also dis¬ 
play pseudogap behavior due to the precursor thermal fluctuations? What is so 
different about the hole-doped cuprates? 

In this review, we shall delve into discussing some of the merits of both ap¬ 
proaches. However, hopefully by the end of this article, it will become clear to 
the reader that we believe that the spin liquid approach provides a more coherent 


a Though clearly important for describing some of the phenomenology, we shall ignore the effects 
arising due to the presence of quenched disorder , which in some cases forbids a true symmetry 
breaking in two spatial dimensions. 

b By ‘topological order’, we mean states with emergent gauge excitations; for states with an energy 
gap, there are non-trivial ground state degeneracies on a torus. 
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picture of the hole-doped side of the phase diagram. 

Let us now outline the scope of this review. There are already a number of 
excellent reviews that discuss much of the previously known phenomenology of the 
pseudogap phase 5,11,16 . Our emphasis here will be on the more recent experiments, 
particularly with regard to the discovery of charge order in the non-Lanthanum- 
based cuprates, and the implications of these experiments on the nature of the high 
temperature pseudogap phase. We shall not explicitly consider experimental reports 
of nematic order 4 : this order can be an ancillary consequence of our low and high 
temperature models of the pseudogap, but we do not believe it is crucial to the basic 
phenomenology. Also, we shall not discuss any of the signatures of time-reversal 
symmetry-breaking that have possibly been seen via polarized neutron scattering 17 
and Kerr-rotation 18 experiments, but surprisingly not in /iSR or NMR experiments. 
While these are clearly interesting experiments that need a theoretical explanation 
in the future, we shall choose not to discuss them any further in this review. 

With the above goals in mind, let us raise a few sharp questions with regard to 
the pseudogap phase, and the charge density wave in particular, that we hope to 
address in the rest of this review: 

• What is the normal state out of which the charge density wave emerges? 
In other words, can the charge density wave help us identify the nature of 
the parent pseudogap state? 

• What is the modulation wavevector, and, the internal symmetry (‘form- 
factor’) of the charge density wave? 

• What is the nature of the onset, and in particular, the temperature depen¬ 
dence of the charge density wave correlations? 

• Does the charge density wave lead to Fermi surface reconstruction that is 
consistent with quantum oscillations at high fields, and if so, which under¬ 
lying Fermi surface does it reconstruct? 

• Is there an underlying quantum-critical point (QCP) that controls the 
physics of the strange-metal region, and if so, what are the two phases 
on either side of this QCP? 

The reader will probably agree that most of these questions, while easy to state, 
might not have a complete answer at present. However, we shall argue that based 
on a body of work done over the past few years, it is at least partially possible 
to answer some of them. Some of the above questions will necessary require us to 
speculate, with the hope that future work will help us address these issues further. 

The rest of this review is organized as follows: In Section 2, we summarize some 
of the recent path-breaking experiments that have helped elucidate the nature of 
the charge-order in the underdoped cuprates. We also define a convenient way of 
expressing the modulation of the charge density in terms of a bond-observable. In 
Section 3, we shall describe the phenomenology associated with the low temperature 
pseudogap in terms of fluctuating orders at zero field and long-range charge-order 
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at high-field, without reference to any particular microscopic theory. In Section 4, 
we shall depart from this ‘classical’ perspective and consider a theory for a ‘frac¬ 
tionalized’ Fermi-liquid (FL*) in the high temperature pseudogap. This will involve 
two alternative routes towards arriving at a description of doped carriers moving 
in a background spin-liquid. In Section 5, we shall attempt to connect the low and 
high temperature descriptions by first reviewing the instabilities of a Fermi-liquid 
interacting via antiferromagnetic (AFM) interactions. This approach will lead to 
results that are only partly consistent with the experiments, but will give crucial 
insight into the structure of the problem. However, it will serve as a nice point of 
departure into the FL* perspective: investigating the weak-coupling instabilities of 
the FL* will lead to good agreement with the current experimental observations. 
We shall conclude with an outlook in the final section. 


2. Review of some recent experiments 

In the past few years, a number of remarkable experiments have given both di¬ 
rect and indirect evidence for the appearance of charge-order in the underdoped 
cuprates. Charge (and spin)-stripes have been known to exist in the Lanthanum- 
based cuprates 11,19 for a long time now. In these materials, close to a hole-doping of 
p= 1/8, the spin and charge stripes are known to have the most pronounced signa¬ 
tures; this is also where the superconducting T c goes all the way to zero. However, 
in this review we shall focus only on the non-La-based cuprates, where charge-order 
is not accompanied by any spin-order. There are important differences between the 
properties of these two ‘classes’ of cuprates. 

A preliminary indication of charge-order in Bi2Sr2CaCu208+(5 came from STM 
experiments, which imaged the periodic modulation in the density of states near the 
vortex cores in the Meissner state 20 and in the high-temperature pseudogap state 21 . 
In YBa 2 Cu 30 y , NMR measurements have been able to detect field-induced long- 
range static charge order, predominantly on the oxygen sites, without any sign of 
spin-order 22,23 . Ultrasound measurements at high-fields 24 have further corrobo¬ 
rated a transition to a long-range charge ordered state. The most clear signatures 
of any charge density wave correlations should come from X-ray. Indeed, incom¬ 
mensurate charge order has now been detected with resonant 25,26 and hard 27 X-ray 
scattering, both at zero and high fields. The wave vector has consistently been found 
to be directed along the copper-oxygen bonds and decreases with increasing hole¬ 
doping 28 , a trend that is the exact opposite of what has been known to be the case 
in the La-based-cuprates 14 . Moreover, there are strong indications that the charge 
modulation is highest on the bonds connecting the Cu sites 29,30 . A conclusion that 
is safe to draw from all of these experiments is that there exist incommensurate 
charge density wave correlations in the CuCU plane, which become static and (rea¬ 
sonably) long-ranged upon the application of a high magnetic field, and in general 
compete with superconductivity. 

Some of the key experimental data that have shed light on the nature of the 
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charge-order in the underdoped cuprates, is reproduced in Fig. 2. Fig. 2a from Ref. 
25 shows a peak in the total intensity at an incommensurate wavevector ps 27t( 0.3, 0), 
in the superconducting phase. In order to study the interplay between charge-order 
and superconductivity, Chang et a/. 27 , studied the onset of scattering intensity as 
a function of temperature both at zero and finite fields (Fig. 2b). They observed a 
gradual onset of charge density wave correlations at Tqdw (where T c < Tqbw < T *; 
Tqdw typically correlates with T**, as introduced earlier), with a subsequent de¬ 
crease in the intensity below T c at zero-field. The feature for T < T c highlights 
the competition between the two order-parameters, which goes away when super¬ 
conductivity is suppressed by a reasonably large magnetic field. The field-induced 
transition in the charge-order has also been investigated using NQR experiments 23 
(Fig. 2c). At moderate fields, such as those used in the X-ray and NQR experiments, 
the charge-order likely develops in the large ‘halos’ surrounding the superconduct¬ 
ing vortices 0 in the mixed-state. The field-induced transition occurs when the halos 
start to overlap, ‘locking’ the orientation of the charge density wave over long dis¬ 
tances. 

Finally, STM measurements show that the primary modulation resides on the 
Cu-0 bonds and has a non-trivial pattern in real space (Fig. 2d,e). The exact 
symmetry properties associated with this modulation have only been unveiled re¬ 
cently 32 , as we shall discuss shortly. 

Let us first prescribe a systematic approach for expressing the charge-modulation 
on a two-dimensional square lattice. The charge density wave, or more appropriately 
bond-density wave (BDW), can be expressed as a bond-observable defined on sites 
L j, 


P H = (4a C i«> = X X 

Q ^ k 


0 ik-(ri-rj 


0 iQiri+rj)/2 


(i) 


where cj a creates an electron at site i with spin tr(=t?i)- The action of various 
symmetry operations become quite transparent for the above parametrization in 
momentum space in terms of a relative momentum, k , and a center-of-mass mo¬ 
mentum Q, as emphasized in Ref. 33. In momentum space, we can write Eq. (1) 
as, 

( C fc-Q/2,a C fc+Q/2,a) = ( 2 ) 

In contrast, the early work of Nayak 34 , and numerous analyses since 35-40 , have used 
the parameterization 

( C fc,cTfc+Q,a) = F Q {k). (3) 

Of course, we can easily go back and forth between (2) and (3). However when 
we write Pg(fe) as a periodic function of k alone, then the equivalent Fg(fe) also 


c In the extreme type-II limit (/c = A/£sc ~^ °°\ A being the London penetration depth and £sc 
the superconducting coherence length), the vortices really correspond to ‘superfluid’ vortices 31 . 
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Fig. 2. (a) Dependence of the CDW intensity at 15 K, showing a finite correlation length (as 

inferred from the width of the peak). Adapted from Ref. 25. (b) Temperature dependence of 
the peak intensity as a function of magnetic field. The onset shows a gradual non-mean-field like 
behavior. At zero-field, the intensity decreases below the superconducting T c , showing competition 
between the two order-parameters. Adapted from Ref. 27. (c) Quadrupole part of the splitting 
of the O line showing a field-induced transition. Adapted from Ref. 23. (d), (e) R-Maps taken at 
150 mV. Adapted from Ref. 29. 


depends upon Q. In particular the state described by Pg(fe) ^ cos(k x ) — cos (fey), 
proposed in Refs. 41,42, is a d-form factor bond density wave which preserves time- 
reversal and which will play an important experimental role shortly. In contrast, 
the state Pg(fe) ^ cos (k x ) — cos (fey), considered by others, is a different state; 
for general Q, it is a mixture of components that are both even and odd under 
time-reversal, and so it is not a useful starting point for a symmetry analysis. 

Conventional charge density waves with Pg(fe) independent of k lead only to 
an on-site charge modulation, with the overall modulation period set by 2tt/\Q\. 
However, this is not the case in the context of the cuprates. The experiments on at 
least two different families of the cuprates (BSSCO and Na-CCOC), which involves 
phase-sensitive STM 32 and X-ray 43 measurements, have now unveiled the form 
factor Pq(k) to be predominantly of a d-wave nature, i.e. Pg(fe) ~ (cos k x — cos k y ). 
In addition, as already mentioned above, almost all the experiments point towards 
a strong evidence for the wavevector Q to be along the Cu-0 bonds, i.e. (±Qo,0) 
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and (0, ±Qo)? where Q o decreases with increasing hole-doping (and is therefore, 
incommensurate). For an underlying ‘large’ putative Fermi surface, which is absent 
in this regime (as known from photoemission experiments, which only see Fermi- 
‘arcs’), these wavevectors would nest regions in the vicinity of, but away from the 
antinodes: ( 71 , 0) and (0,7r). 

On the other hand, there is a recent X-ray observation 44 in the La-based 
cuprates, which seems to indicate a dominant s'—form factor (where Pg(fe) ^ 
(cos k x -fcos fc y )), and this has been ascribed to the presence of magnetic stripe order 
in these compounds 45,46 . The nature of the bond-density waves seen in the two dif¬ 
ferent ‘classes’ of cuprate superconductors is therefore qualitatively different—both 
in terms of the predominant form-factor, and, in terms of the doping dependence 
of the wave vector. 

In Fig. 3 we provide an illustration of unidirectional BDWs with different com¬ 
ponents of Pg(fe) in real space for both commensurate as well as incommensurate 
wavevectors. It is not a coincidence that the patterns of Fig. 3 (c), (f) and the data 
of Kohsaka et al. 29 in Fig. 2(d), (e) look remarkably similar. This comparison has 
recently been carried out to a remarkable degree of precision in the phase-sensitive 
work of Fujita et al. 32 , which has pinned down the symmetry to be predominantly 
d-wave. 

A key question now is whether the presence of the charge density wave order may 
also explain the quantum oscillations observed at low temperatures and high mag¬ 
netic fields 48,49 . We shall discuss this topic at some length in the next section, but 
it is worthwhile to also review some of the most striking experimental aspects here. 
The first observation of quantum-oscillations 48 provided clear signatures of a recon¬ 
structed pocket with an area approximately 2% of the Brillouin zone (corresponding 
to a frequency ^ 530T; see Fig. 4a). This is fundamentally different from the ob¬ 
servations in extremely overdoped cuprates, where quantum-oscillations observed a 
‘large’ Fermi surface 50 . The key question then was to figure out the density-wave 
order responsible for giving rise to the reconstruction in the underdoped cuprates. 
An interesting fact that was realized soon after the initial discovery of oscillations, 
is related to the nature of carriers in the system. In spite of doping holes into the 
system, transport-properties were most consistent with the presence of primarily 
electron-like quasiparticles, as evidenced by the negative Hall coefficient 51 (Fig. 4b). 
This prompted the search for a small electron-pocket , that would explain both of 
these observations. 

More recently, by extracting the effective cyclotron mass of the electron-like 
quasiparticles d and by noting the critical dopings where they appear to diverge 
(Fig. 4c), the location of the underlying quantum critical points, corresponding to the 
onset of charge-order, have been determined 49,52 . These QCPs also coincide with 
the dopings around which the T C (B)- domes, as deduced from the resistive transition 


d By studying the temperature dependence of the oscillation amplitudes and fitting to the Lifshitz- 
Kosevich form. 
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(d) (e) (f) 


Fig. 3. Real-space visualization of unidirectional BDW (charge-stripe) with components: P s (k) = 
P s ,P s /(k) = P s /(cosk x + cos ky), and P d {k) = P d (cosk x — cosk y ). (a)-(c) we plot the charge 
modulation for a commensurate wavevector, Q = 27t(|, 0), while (d)-(f) plot the same quantity 
for an incommensurate wavevector, Q = 27t( 0.3,0). The parameters used are: (a), (d) P s = 
0, P s / =1 ,P d = 0. (b), (e) P s = 1, P s t = l,P d = 0. (c),(f) P s = 0, P s / =0,P d = 1. We have 
included phases in the definitions of P s s / d in order to make the charge distribution bond-centered 
for the case of the commensurate wavevector; however, other choices of the phases are also allowed. 
Adapted from Ref. 47. 

in a field (£>), are centered about. We direct the readers to some recent review 
articles for further details regarding the quantum oscillation experiments 49,53 . 

(a) (b) (c) 





Fig. 4. (a) The oscillatory part of the Hall resistance (after subtracting the monotonic background 

at T = 2 K), as a function of inverse magnetic field, 1/B. Adapted from Ref. 48. (b) The Hall 
coefficient as a function of B at different temperatures. Adapted from Ref. 51. (c) The blue curves 
represent T c , at different B and are supposed to highlight the connection between effective mass 
enhancement and the strength of superconductivity. Adapted from Ref. 52. 
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3. Phenomenology at lower temperatures 

After the above review of some of the key experiments, let us start addressing the 
questions that we had raised at the end of the introduction. Here we shall present 
semi-phenomenological models of the ordering instabilities at low temperatures. 

3.1. Onset of charge-order and superconductivity 

As mentioned in the previous section, a peculiar feature associated with charge- 
order in the underdoped cuprates, is the highly non mean-field-like nature of the 
onset, at a temperature that is typically lower than the pseudogap temperature, T*. 
In fact this behavior is reminiscent of the onset of antiferromagnetism in the parent 
insulating antiferromagnetic state of La 2 Cu 04 54 . We now understand this gradual 
onset of intensity 6 in the structure factor over a broad range of temperatures in 
terms of the angular fluctuations of an 0(3) field (corresponding to the Neel order) 
in two spatial dimensions 55 . 

The above observation then begets the question if the onset of charge-order 
can be understood in terms of the angular-fluctuations of an order-parameter 56 , 
and if so, can one write down an effective field theory that describes these fluc¬ 
tuations. Let us first note the symmetries associated with the above problem: 
0(2) xO(2)xO(2)XZ 2 . These correspond to, charge-conservation, x— translations, 
y —translations and the x y symmetries. In addition, we use the fact that the 
system also preserves time-reversal and inversion symmetries. 


Pseudogap 

¥ 


Superconductor 

v 


T <T C 



Fig. 5. The angular fluctuations of the six-component order-parameter, describe the non-mean 
field onset of charge density wave correlations below Tcdw- Due to a finite g > 0, fluctuates only 
along the superconducting directions for T < T c , the Kosterlitz-Thouless transition temperature 
corresponding to the 0(2) component associated with T. For w < 0, the Z 2 invariance associated 
with x <ry can be broken spontaneously. 


The onset can then be best understood in terms of the angular fluctuations 


The finite transition temperature arises solely due to the non-zero interlayer coupling. 
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associated with a six-component order parameter 57 : 

ft = (ReT, ImT, Re^, Im$ x , Re<h y , Re<h y ), 


( 4 ) 


where T represents the superconducting component and <& x , y represent the incom¬ 
mensurate charge density wave order parameters along x and y directions^ The 
fluctuations are assumed to be of a purely classical and thermal nature and are 
postulated to be described by the following non-linear cr-model (NLcrM), 


S c \ 

~T 


Z = J vn(r)5(rt 2 - 1 ) exp ^ - 
5 c i = y J d 2 r |VM /| 2 + A(| V $ K | 2 + | V$ y | 2 ) 


+ 5 (|^| 2 + |^| 2 ) + U ;( I ^| 4 + I ^| 4 ) 


( 5 ) 

(6) 


In the above, ps and control the helicity moduli of the superconducting and of 
the density wave order respectively. The coupling g breaks the symmetry explicitly 
between the T and <$> y directions. It sets the relative energetic cost of ordering 
between the superconducting and density wave directions. We take g > 0, so 
that superconductivity is preferred at low temperatures (figure 5). The coupling w 
imposes the square lattice point group symmetry on the density wave order. 

Note that within this framework, in the intermediate temperature range (T < 
Tcdw) over which the charge order correlations have been observed, we can choose 
not to include the effect of antiferromagnetic correlations explicitly; they can be 
absorbed into the phenomenological parameters of the above classical model. How¬ 
ever, for the onset of the pseudogap at T*, these can’t be ignored. 

It is possible to add higher-order terms to the above action, without significant 
qualitative differences in the results. The above action was studied using (classical) 
Monte-Carlo simulations and within a ‘1 /TV’ expansion (where N corresponds to 
the number of components of fi; N = 6 for the theory in Eq.6) 57 . The results for 
the structure factor and comparison to experiments on YBCO are shown in Fig. 
6. The agreement with experimental data is very good, except at low and high 
temperatures. In the absence of impurities, that would otherwise pin the charge- 
order and give rise to a finite intensity at low temperatures, the structure-factor 
goes to zero. On the other hand, at higher temperatures, the onset of amplitude 
fluctuations, which are ignored in the above setup, would also decrease the intensity. 

The theory in Eq. (6) has also been applied to measurements of the diamagnetic 
susceptibility 58,59 in YBCO. This measures the superconducting T component of 
the 0(6) order parameter, complementing the X-ray measurements of the charge 
order represented by <& x , <& y . It was found 60 that the same set of parameters 
can provide a reasonable simultaneous fit to both the X-ray and diamagnetism 


f Note that the region in the vicinity of $7 = 0 is forbidden, i.e. the amplitude fluctuations of |fi| 
are completely forbidden below Tcdw and probably onset at an even higher temperature. 
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Fig. 6. Comparison of X-ray data for the CDW scattering intensity to the structure-factor cal¬ 
culated in Monte Carlo simulations of Z for different parameter sets. Adapted from Ref. 57. 


measurements. Thus the 0(6) theory provides a surprising quantitative connection 
between two very different classes of experiments, and this supports its validity as 
a theory of fluctuating orders in the low temperature pseudogap regime. 


3.2. Quantum oscillations 

Some of the most interesting aspects of the underdoped cuprates involve the na¬ 
ture of the fermionic excitations. However, in the above description so far, we 
have ignored the Fermionic degrees of freedom altogether. In order to study these 
questions, here we couple the underlying electrons to d-wave superconducting and 
d-form-factor charge density wave order parameters and explore the fermionic ex¬ 
citations at low temperatures and high magnetic fields. As discussed in Section 2, 
this is the regime where quantum oscillations have been observed. It is not entirely 
clear what the upper critical field, H c2 , is for the underdoped cuprates and if the 
fields, B ~ 40T, where oscillations are first observed are already higher than id c2 . 
Moreover, whether the application of large fields completely suppresses the pseudo¬ 
gap is also a topic of ongoing debate. These are questions that need to be resolved 
in order to realize what gets reconstructed to give rise to the small electron-pocket. 

It is nevertheless a reasonable starting point to study the problem of reconstruc¬ 
tion on top of a ‘large’ Fermi surface in the presence of a (sufficiently) long-ranged 
BDW S . 

The spectral function, in the unfolded Brillouin zone, for a large Fermi surface 
in the presence of long-range h BDW is shown in Fig. 7a. A small nodal electron- 


g Such that the cyclotron-radius associated with an applied magnetic field is smaller than the 
correlation length of the BDW. 

h Such long-ranged BDW is likely to be present only in a strong magnetic field and will not be seen 
in photoemission experiments. 
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(a) (b) 1 



Fig. 7. (a) The electron spectral function in the presence of long-range bidirectional BDW (with 

wavevectors Q ± and Q 2 ) at zero magnetic field (in the unfolded Brillouin zone). In black, the 
Fermi surface used for our computation. The BDW causes reconstruction of the Fermi surface 
and the formation of an electron-like pocket (red) and two hole-like pockets (blue), (b) Electron 
spectral function in the presence of fluctuating superconducting and BDW correlations. Adapted 
from Ref. 61. 


pocket 62 can be seen (shown in red) 61 , which is likely responsible for giving rise to 
the oscillations corresponding to ^ 530T. A significant prediction of our computa¬ 
tion is that an invariable consequence of reconstructing a large Fermi surface is the 
presence of small /io/e-pockets closer to the antinodes (shown in blue). There have 
been recent reports for the existence of these hole-pockets 63 at frequencies close to 
those predicted by our analysis. 61 

We have also extended the above analysis of reconstruction of the large Fermi 
surface to the case of fluctuating order: this would apply at zero field in the same 
lower temperature regime where we described X-ray and diamagnetism experiments 
in Section 3.1. We describe the fluctuating order by the 0(6) model, and compute 
the fermion self energy at the one loop level. Typical results appear in Fig. 7b, 
showing ‘Fermi arc’ spectra in the nodal regions. In photoemission experiments, 
such Fermi arc spectra are also seen at significantly higher temperatures, where 
the 0(6) model of fluctuating order cannot be reasonably applied. We believe a 
rather different model is necessary at higher temperatures, and we will turn to this 
in Section. 4. 

Recent experiments have also detected a divergence in the effective cyclotron 
mass, ra*, as extracted from quantum-oscillations, near optimal doping 52 . It has 
been argued 64 that the critical like divergence arises from the angle-dependent 
Fermi-velocity, vp(0 ), along the nodal electron-pocket, and is in fact dominated by 
the corners of the pocket 1 . However, signatures of this divergence are absent in 


^ince m* = T- f . 

2-7T J Vp{0) 
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any measurement performed at zero field, that is dominated by the ‘light’ nodal 
quasiparticles, such as transport. 

Despite the agreement between the computation of the large Fermi surface re¬ 
construction by the BDW and the quantum oscillation experiments, there are a 
number of outstanding puzzles that remain. Notice that in addition to the closed 
pockets in Fig. 7(a), there are a number of open Fermi-sheets, that are bound to be 
present as long as the BDW potential is small. Even though these don’t contribute 
to oscillations, the large density of states along these sheets would contribute to var¬ 
ious thermodynamic measurements. However, recent measurements of the Knight- 
shift 65 and specific heat 66 in fields of upto ~ 45 T do not seem to be consistent 
with the availability of such large density of states at the Fermi-level. One could 
argue that if the BDW potential were to be arbitrarily large, it could gap out these 
portions of the Fermi surface as well. However, in the presence of such strong in¬ 
commensurate density-wave order, quantum oscillations are likely to be disrupted 
completely 67 . 

The above thermodynamics experiments seem to suggest that whatever is caus¬ 
ing the pseudogap to open up close to the antinodes remains unaffected at low to 
intermediate fields. This leaves us with two options: (i) The reconstruction oc¬ 
curs on top of a ‘pseudogapped’ Fermi surface, so that it is really the arc that 
gets reconstructed. We shall revisit this scenario further in Section 5.3, or, (ii) 
The reconstruction happens in a ‘vortex-liquid’ phase, where the vortices associ¬ 
ated with d -wave superconductivity are fluctuating 68 . However, within either of 
these two scenarios, the nodal electron-pocket is likely to survive, since this region 
of the Fermi surface remains unaffected, both by d-wave superconductivity and by 
the pseudogap. 

Another issue that remains to be resolved with regard to the nature of the 
charge-order responsible for the oscillations at high fields is as follows. Experiments 
at zero-fields seem to suggest the presence of a stripe-like order, while the nodal 
electron-pocket scenario requires a nearly checkerboard order. In order to relate 
these two regimes, there could be two possibilities: ( i ) There is a field induced 
transition from a stripe-like BDW at zero, or, low-fields to a nearly checkerboard- 
like BDW at high-fields, or, (ii) The system actually consists of differently oriented 
domains with stripe-like order, either within each Cu 02 plane, or within different 
planes (but with a strong interlayer tunneling) 69 , where the typical domain sizes, 
^domain, are smaller than the cyclotron-radius. Future experiments will hopefully 
be able to shed light on this issue. 


4. Fractionalized Fermi liquid at higher temperatures 

We now turn to the higher temperature regime just below T*, significantly above 
temperatures at which there is any indication of appreciable fluctuating order, and 
so the analysis of the previous section cannot directly apply. Here, as we have 
discussed earlier, the primary experimental indications of the pseudogap are the 
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suppression in the spin susceptibility, and the appearance of ‘Fermi arc’ spectra in 
the electron spectral function. 

We review a model of this pseudogap metal as a ‘fractionalized Fermi liquid’ 
(FL*). While related states have been derived by many different theoretical meth¬ 
ods 6,70-72 , we begin with a simple toy model description 73 which generalizes the 
quantum dimer models 74,75 of insulating spin liquids. (We believe that recent re¬ 
sults from dynamical mean field theory 76,77 also point to a low energy effective 
theory of the pseudogap in terms of such a dimer model 73 ). Consider a doped an- 
tiferromagnet with p holes per unit cell, as shown in Fig. 8a. Note that there are 




<E3H< ■ 





CD <3 E3> 


(a) Spin liquid with holons 


Fig. 8. (a) A density of p holes per unit cell in an ordered antiferromagnet. The holes are the 

open circles, while the arrows denote the orientation of the electron spins on the remaining sites. 
Note that there is a density of 1 + p holes with respect to the band insulator in which there are 
2 electrons on each site, (b) Spin liquid with a density p of spinless holons of charge +e. The 
ellipses represent spin-singlet pairs of electrons in the state (|N) — |4T)) /v^2- 


p holes with respect to the antiferromagnet, but the number of holes with respect 
to the filled band insulator is 1 + p. As p increases, we imagine that the antiferro¬ 
magnetic order is quickly destroyed, and we obtain a state with p holes in a spin 
liquid background as shown in Fig. 8b. Now all the spins of the antiferromagnet 
have been paired into singlets, and these singlets can resonate with each other 78 . 
The holes in this spin liquid background carry no spin and charge +e as they move 
around: so these are ‘holons’. This doped spin liquid state also has neutral 5 = 1/2 
excitations known as ‘spinons’, as shown in Fig. 9a. 

So far we have just described the well-known structure of a doped spin liquid. To 
obtain the FL*, we need one further step: suppose there is an attractive potential 
between the holons and spinons so that they form bound states with charge +e 
and spin S = 1/2 (a rationale for this attractive potential is given in the caption of 
Fig. 9a). For simplicity, let us take the strong coupling limit in which this bound 
state involves only nearest neighbor sites, and so can be represented by a dimer, as 
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Fig. 9. (a) An excited state of the doped spin liquid, with neutral S = 1/2 spinon excitations 

represented by the arrows. Both the holons and the spinons can move independently in the spin 
liquid background. Note the holon and spinon on nearest neighbors in the middle of the figure: 
the electron can easily exchange positions between these sites, and this can lead to a holon-spinon 
bound state in which the electron is in a bonding orbital between the two sites, (b) The FL* state: 
All holons have paid the energy cost to create a spinon from the spin liquid, and the resulting 
bound state is represented by the green dimers. The green dimers are fermions carrying charge 
+e and spin S = 1/2, and their motion in the spin liquid background realizes a metal with Fermi 
surfaces enclosing an area equivalent to a density of p fermions of spin S = 1/2. 


shown-i in Fig. 9b. If the binding energy is large enough, every holon will pay the 
energy cost needed to create a spinon out of the spin liquid background, and we 
obtain a spin liquid doped only by fermionic dimers carrying charge +e and spin 
S'= 1/2. 

The dynamics of this FL* is now described by a quantum dimer model 73 in 
which the fermionic dimers and the original spin singlet pairs resonate with each 
other. This leads to motion of the fermionic dimers, and as there is a dilute gas 
of them, we can expect them to form a metallic state with a Fermi surface. The 
quasiparticles near the Fermi surface will now have the same quantum numbers as 
in a Fermi liquid, but they differ in key aspects: 

• The volume enclosed by the Fermi surface in a FL* is equivalent to a density 
of p spin 5 = 1/2 particles of charge +e. In contrast, in a Landau Fermi 
liquid the equivalent density would be 1 + p. 

• The FL* has a second independent sector of low energy ‘topological’ exci¬ 
tations. In the quantum dimer model, these are associated with the reso¬ 
nances of the spin-singlet dimers, which are represented in analytical for¬ 
mulations by emergent gauge fields. 

• The quasiparticles near the Fermi surface are neutral under the emergent 


J The bound state is represented in Fig. 9b by an arrow, and this represents an electron which is 
in a bonding orbital between the two sites. In a three-band model, this can be associated with 
an electron on the O orbital between the Cu orbitals. However, our analysis does not require the 
three-band model, and applies equally to the one-band Hubbard model. 
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gauge field. In the quantum dimer model presented here, this is a conse¬ 
quence of the fact that the fermions are themselves also dimers. Indeed, 
the fermions couple to the gauge field non-minimally as dipoles. This cou¬ 
pling is not strong enough to disrupt the quasiparticle nature of the Fermi 
surface excitations. 

The dimer model picture demonstrates the inevitability of the above characteristics 
of the FL*: the fermionic dimers have density p, while the remaining sites are 
required to have spin-singlet dimers representing the gauge field. In Ref. 79,80, 
general arguments have been given which demonstrate that any violation by a Fermi 
surface of the 1 +p Luttinger volume requires the presence of low energy topological 
excitations. 

The following subsections will move beyond this dimer toy-model picture of 
the FL* phase, and review field theoretic analyses which yield a metallic state 
with the same basic characteristics. Such analyses allow us to compute the Fermi 
surface structure, understand potential low T instabilities of the FL* metal, and 
also connect to other phases in the global cuprate phase diagram. 


4.1. Field-theoretic description 


We begin our description of the FL* metal by first reviewing a standard description 
for the evolution of antiferromagnetism in metallic phases of the one-band Hub¬ 
bard model. This model is widely believed to contain the essential physics of the 
cuprates 78 , for electrons hopping on the sites of a square lattice, 




hubbard 


= H t + Hr 


i<j 


Hu =UJ2 ("it ~ \) ( «4 - 



(7) 

.t c . 

(8) 

1\ 


' 2 ) 5 

(9) 


where Uj represent the hopping parameters, U is the on-site Coulomb repul¬ 
sion, (i represents the chemical potential and a =t,l are the spin-indices. The 
Fermions satisfy the standard anti-commutation relations, {Q a ,c^} = SijS a /s and 
{q a ,c^} = 0. This completely defines the problem at hand, that has eluded an 
exact solution. 

Anticipating the metallic state to be in the vicinity of an antiferromagnetic 
instability close to half-filling, we use the exact operator equation, 


U 





U 

7 ’ 


( 10 ) 


valid on each site, i, and where Si = c\ a cr a ^Cip/2. Upon decoupling the interaction 
via a Hubbard-Stratanovich transformation, we obtain 


exp 




J VJ z (r)ex - Yi 


dr 
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We can now integrate out the fermions from the action and look for the saddle 
point of the resulting action for J This leads to the Neel state with a wavevector, 
K = (7r,7r). In the long-wavelength limit, it is then useful to introduce a field, 

Ji = <Ae* K " ri . (12) 


This is then the familiar route towards arriving at the “spin-fermion” model, 


Z = 

s = 

+ 


J Vc a V(f exp (-S), 

dT C L( 5 t - £k)Cka 4* Vi ' 

J L i„ 

/ 


Oif3Ci(3£ 


;k -n 


k 

dr d 2 r 


— (V r (y?) 2 + ~{d T (f) 2 + -ifi 2 + ~^( 1 P 2 ) 2 


(13) 


(14) 


The first term in S is just the kinetic energy of the fermions, and the second term is 


(0)¥=o 

0 



0 


AF Metal with 
hole pockets 


0 



AF Metal with 
electron and 
hole pockets 


{<f) = 0 



Metal with 
“large” Fermi 
surface 


5 


Fig. 10. The evolution from a ‘large’ Fermi surface (for (cp) = 0, s > s c ) to a reconstructed Fermi 
surface with electron (red) and hole (blue) pockets in the presence of antiferromagnetic (AF) order 
with a finite ((f) / 0 (s < s c ) across a SDW QCP at s = s c . By moving sufficiently deep inside 
the SDW phase, one has a metallic state with only hole-pockets in the presence of a large (cp). 


the “Yukawa” coupling term that leads to the spin-fluctuation mediated scattering 
of fermions from one portion of the Fermi surface to the other. The momentum of 
the boson, <p, is small, and v represents a characteristic spin-wave velocity. The last 
line in the above action is just the cp 4 — field theory for aniV = 3 order-parameter. 

At mean-field level, as we tune the value of 8, we go from a metallic Fermi-liquid 
phase with a ‘large’ Fermi surface (s > s c ) to a metal with reconstructed electron 
and hole pockets (s < s c ). The gap due to the SDW order parameter opens up at 
the ‘hot-spots’, where = £k+K = 0. With a decreasing 8, as (cp) increases in 
magnitude, only the hole pockets remain (Fig. 10). 
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So far, apart from possible exotic behavior at quantum critical points, we have 
only obtained conventional Fermi liquid phases. To obtain the FL*, we have to 
switch from the above description of SDW order by a “soft-spin” cp field with large 
amplitude fluctuations, to a “hard-spin” perspective in which we have primarily 
angular fluctuations of the antiferromagnetic order. So we replace ip by a unit 
vector field n, 

(p => n , n 2 = 1. (15) 

The key utility of such a formulation in terms of the unit-length n field is that it 
allows us to consider a route to “quantum-disordering” the SDW order in which 
topological “hedgehog” tunneling events are suppressed 81 ; in contrast, the pertur¬ 
bative analysis of the soft-spin theory necessarily proliferates hedgehogs at the zeros 
of the (p field. In such a hard-spin theory, we argue that the evolution of phases in 
Fig. 10 can be replaced by the more exotic route shown in Fig. 11. Now there is an 



g 


Fig. 11. Alternative route to Fig. 10 for the loss of antiferromagnetic order in a metal in a theory 
with angular fluctuations of the fixed-length order n. The middle phase is FL*. The evolution 
from FL* to the large Fermi surface Fermi liquid on the right is addressed in Ref. 82. 


intermediate phase in which there are hole pockets enclosing a volume equivalent 
to fermions with density p, but without antiferromagnetic order: this is the FL* 
phase. 

Let us turn to an explicit presentation of the theory of a metal with angular fluc¬ 
tuations of antiferromagnetic order. We are only interested in the long-wavelength 
fluctuations of rq while retaining the full lattice dispersions for the fermions. Ap¬ 
plying (15) to (14) we obtain the imaginary time Lagrangian, C = Cf + C n + £/ n , 
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with 


c f = X/ 


{d T t' 


Cj ex H - h. c. 


£„ = 2g / d2r +v2 ( Vn ) 2 ]’ 

C fn = Xj2e lKr ' rij • 4 a a a0 c i0 . 


(16) 

(17) 

(18) 


Now s has been replaced by g to tune the strength of quantum fluctuations associ¬ 
ated with the AFM order parameter. We assume g is a generic coupling measuring 
the degree of frustration in the insulating antiferromagnet, which can drive the 
insulator into a non-magnetic state with valence bond solid (VBS) order across a 
deconfined quantum critical point 83 . It is therefore useful to discuss the phase di¬ 
agram of a general class of frustrated doped antiferromagnets, as a function of the 
coupling g and the charge density p, as shown in Fig. 12. 



Fig. 12. Adapted from Ref. 47. The vertical axis represents some frustrating coupling in an 
insulating antiferromagnet which can drive it into a valence bond solid (VBS) across a deconfined 
quantum critical point at g = g c . 


For g < g c , the above model has long-range antiferromagnetic order, (n) ^ 0 
(with a correlation length, £ —>■ oo), which reconstructs the large Fermi surface. The 
interesting regime is when g > g c , which we argue is relevant to the physics of the 
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non-La-based cuprates. The deconfined critical point g = g c and pm 0 is expressed 
not in terms of the n fields, but in terms of the spinor z a , with rq = z* a a a ^Zi^ and 
a U(l) gauge field A^. At the same time 72 , one tranforms the underlying electrons, 
to a new set of spinless fermions, V’zp, 


(19) 


Co = R a P ^ipi where 



( 20 ) 


is a spacetime dependent SU(2) matrix (52 a \zi a \ 2 = 1) and the fermions i/ji p carry 
opposite charges p = ±1 under the same emergent U(l) gauge transformation. Note 
that the above parametrization introduces a large SU(2) ‘gauge’ redundancy and it 
is possible to embed the U(l) gauge-theory within the larger SU(2) gauge structure. 
This has been used to describe the QCP around optimal doping recently, in terms 
of a ‘Higgs’-transition without any broken symmetries 82,84 

We begin our discussion of the phase-diagram in Fig. 12 by first considering the 
effect of non-zero p at the critical coupling g = g c . We are not interested here in 
the very small values of p at which the ^p fermions get effectively localized 71 and 
exclude this regime in the phase-diagram (shown as the grey-shaded region). At 
higher hole densities, the holons form a Fermi surface, and this can then quench the 
A p fluctuations via Landau damping 85,86 ; the holon Fermi surface is then stable k 
against confinement 88,89 and we obtain the holon metal, also referred to as the U(l) 
algebraic charge liquid , shown in Fig. 12. 

Let us now turn to g > g c . Now there is at least one additional length scale, the 
spin correlation length, £. This length should be compared with the spacing between 
the holons ~ 1/y/p. When £ <C 1 /y/p, we have to first consider the influence of a 
non-zero £ on the gauge theory of the deconfined critical point. As described in 
Refs. 83,90, in this regime there is a crossover to a Coulomb phase in which the 
An field mediates a logarithmic Coulomb force. This Coulomb force binds the 
and V^p quanta into gauge-neutral fermions 71,72,87 (an additional attractive force is 
also provided by the “Shraiman-Siggia term” 71,72,9 i). At longer scales, these gauge- 
neutral fermions start to notice each other via the Pauli principle, and so they form 
Fermi surfaces leading to the FL* state of interest here. 

As discussed above, the emergent photon gives rise to binding of the xjj p fermions 
and the spinons into gauge-neutral objects. However, there are two such combi¬ 
nations, 


R ia ^ Zi a 'lpi-\-i Gi a ~ £ a pZip'lpi— , 


( 21 ) 


where e a p is the unit antisymmetric tensor. Note that unlike the bound-state 
Fermionic dimers in the previous discussion of the toy-model, the above bound- 
states arise from spinons and holons on the same site. The physical electronic 


k At low temperatures, the holons can pair to form a composite Boson that is neutral under the 
A p field, condensation of which leads to the holon superconductor 87 . 
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operator has a non-zero overlap with both of these combinations, 


Qa — Z(Fi a H~ Gia'ji 


( 22 ) 


where Z is a quasiparticle renormalization factor that is nonlocal over £. Over 
distances that are larger than £, where there is no net AFM order, the F a and G a 
fermions preferentially, but not exclusively, reside on the different sublattice sites. 

Based on symmetry considerations alone, we can write the following effective 
Hamiltonian 72 for a FL* metal realized by Fi a and Gi a , 

= -E + G h G ja) + A e iK ^(FlF ia - GlGia) 

i,j i 

~'FU F lGja + GlF ja ). (23) 

iJ 


Once again, the represent the hopping matrices corresponding to a large Fermi 
surface, A represents the potential due to the local AFM order (at distances shorter 
than £). The terms proportional to represent the analogs of the “Shraiman- 
Siggia” (SS) terms 91 which couple the F and G particles. When tij ^ 0, one 
obtains hole-like pockets centered away from (7 t/2 ,7 t/2 ). 

The Green’s function for the electronic operator is given by 72 , 


G c (k,u) 


_Z 2 _ 

U ~ ^ 2 /[^ - ik+ k] 


(24) 


where the dispersions, are, 


Cfe" £k T Ski £k £k 5 (25) 

Ek = — 2£i(cos(fca;) + cos(ky)) — At2 cos (k x ) cos (ky) 

—2ts(cos(2k x ) + cos(2 k y )) — /i, (26) 

£k = -to - ii(cos (k x ) + cos(k y )). (27) 


The spectral function corresponding to the above Green’s function, along with 
comparison to experimental photoemission spectra, is shown in the left panel of 
Fig. 13. Note that the theoretical ‘small’ pocket Fermi surfaces have an area given 
by p, and not (1 +p). This gives rise to a violation of Luttinger’s theorem, which is 
however allowed because of the presence of ‘topological-order’ due to the background 
U(l) spin-liquid. However these are pockets of hole-like quasiparticles, which are 
gauge-neutral under the emergent U(l) field, Ay. 

The FL* is a conventional Fermi-liquid phase, as far as its transport, thermo¬ 
dynamic or spectral properties are concerned. However, as noted above, the subtle 
difference arises from the presence of topological order. The spectral weight on 
the ‘back-side’ of these pockets is suppressed strongly. It has been argued that 
the Fermi-arcs in the pseudogap regime could in fact be modeled as pockets whose 
back-sides are (almost) invisible 6,10 . However, clear signatures of the ‘back-sides’ 
are currently lacking in experiments. 
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Fig. 13. On the left is the spectral function of the FL* state computed in Ref. 72. On the right 
are photoemission measurements of the pseudogap phase by Shen et al. in Ref. 92. 

5. Connecting the high and low T regimes 

We will now connect our description of the FL* model in Section 4 to the charge 
order and superconducting instabilities presented in Section 3. 

For theoretical purposes, and also with the aim of reviewing part work, it is useful 
to first present the instabilities of a conventional Fermi liquid state in Section 5.1. 
It turns out that these instabilities lead to some discrepancies with experimental 
observations. We will then argue in Section 5.3 that these discrepancies are naturally 
resolved by replacing the Fermi liquid by the FL* phase as the parent phase of the 
orderings at low T. 




5.1. Ordering instabilities of a Fermi liquid with 
antiferromagnetic correlations 

In this subsection, we return to the soft-spin theory in Eq. (10) for the properties 
of a Landau Fermi liquid in the presence of amplitude fluctuations of a SDW order. 
Starting from Eq. (10) it is useful to focus on the low energy theory for fermions 
near the ‘hot spots’ on the Fermi surface, shown in Fig. 14. Such a Lagrangian is 
given by, 


t-'sdw — H - H - 


Cf = vlr (dr - iv? ■ + 47(9r - iv? ■ V)^ 2 ™ 

£fv = <p ■ (V’i^rp'V’v + 


(28) 


and where C^ part of the action was already written down in Eq.14. In the above, 
^m, v m re p resen t the velocities at the hot-spots labelled ‘A’, ‘C’, or, ‘B’, ’D’ within 
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Fig. 14. A cartoon of a large Fermi surface, with the filled electronic states in blue. The red and 
green circles represent the ‘hot-spots’ connected by K (shown as wavy lines). 


each hot-spot pair (denoted by ‘ra’). 

Before proceeding any further, let us first digress for a bit and discuss an inter¬ 
esting (exact) symmetry associated with the exchange-interactions 5,93,94 , that also 
manifests itself as an emergent symmetry 42 of the hot-spot theory. To start with, 
consider just the nearest neighbor Heisenberg term, 

H.J = Y^J,jS, - S r (29) 

i<j 

where Si is expressed in terms of the underlying electronic operators in the usual 
way. If we, for instance, make the transformation, c\^ -A C 4 , cj^ -a —c^, then 
Sf = — riii ~Sf and similarly for Sf . The above particle-hole transformation 

carried out on every site leaves the Hamiltonian invariant. 

It is possible to generalize the above transformations to arbitrary SU( 2 ) rotations 
in particle-hole space. In order to do this, let us introduce the Nambu spinors, 



Then, Hj can be re-expressed as, 

Hj = l Y, A (*L <r*P*WaJ ■ , (31) 

i<j ' ' ' ' 

where a, b are Nambu indices. It is now a trivial matter to see that Hj is invariant 
under independent SU( 2 ) pseudospin transformations acting on each lattice site of 
the form -A Ui^iab- 

We note that this pseudospin symmetry is broken explicitly by terms involving 
Coulomb repulsion and a finite chemical potential (except when we are at a fine- 
tuned, half-filled state). However, an interesting feature of the low-energy critical 
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theory in Eq.28 is that the SU(2) symmetry re-emerges. To see this, note that in 
Eq.28, if we carry out the following transformations, 


V’lt V’U, V>u “V’lt, ^2t ^24., ->■ -^2t, ¥> -> V, ( 32 ) 


then is invariant because, 


V’ItC Ul ' • V)^u = V’ul'Wi • = ^4,(^1 • V)'0 U , (33) 


where we have used integration-by-parts. The remaining terms can be simplified 
similarly. Note that it is crucial for the curvature-terms to be absent in order for this 
symmetry to manifest itself. The SDW theory with 4 pairs of hot-spots therefore 
has an emergent [SU(2)] 4 symmetry, in addition to the usual spin SU(2) rotation 
symmetry. 

Returning to the problem of a Fermi-liquid in the presence of SDW fluctuations, 
we recall a result that has now been known for almost three decades. The spin- 
fluctuations give rise to a pairing-glue in the d x 2_ y 2 channel 95,96 , as can be seen 
by explicitly evaluating the ‘Cooperon’ diagram (Fig. 15 (a), (b)). The underlying 
reason for the attraction in the sign-changing d-wave channel can be traced back to 
the BCS gap-equation for the pairing gap A*,: 



(34) 


P 


x(p — ft) is the SDW susceptibility, peaked at p — ft = K, the momentum carried 
by the SDW order-parameter and E p = is the dispersion for the Bogoli- 

ubov quasiparticles. In order for the above equation to have a non-trivial solution, 
sgn(A/c) = —sgn(Ap), whenever p — k = K. 

We can now use the SU(2) symmetry to carry out a particle-hole transformation 
on only the electrons in the vicinity of ‘B’ and ‘D’ (Fig. 14). This then suggests 
that the same antiferromagnetic fluctuations would give rise to ‘pairing’ in the 
particle-hole channel as well, as shown by the ‘diffuson’ diagram in Fig. 15(c). The 
particle-hole condensate, (<4 ,_q/ 2 a Cfe + Q/ 2)a ), now carries a finite center-of-mass 
momentum, \Q\ = 2 kp and there is an internal d x 2_ y 2 form-factor, Vd (Fig. 15 (d)). 
Note that this wavevector will always be oriented along a diagonal direction, i.e. it 
is of the form Q = (Q o, ±Qo) where . 

For the hot-spot theory, the superconducting and density-wave instabilities are 
exactly degenerate. However, a finite curvature which breaks the SU(2) symmetry, 
will lift this degeracy. While the superconducting instability remains unaffected 
(since ±k are always ‘nested’), the density-wave instability becomes sub-leading to 
superconductivity. 

The most natural question, in light of our discussion so far, is as follows: Is 
there any relation between the density-wave instability found above in the vicinity 
of a SDW quantum-critical point, to the state unveiled in the underdoped cuprates? 
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(C) 





Fig. 15. (a) ‘Cooperon’-diagram in the particle-particle channel leads to (b) d -wave pairing with 

a relative phase-shift between pairs £ A-B’ and ‘C-D\ (c) ‘Diffuson’-diagram in the particle-hole 
channel, obtained from (a) by carrying out the SU(2) transformation on only half of the hot-spot 
pairs leads to (d) d-form factor density wave with Q = 2k f and a relative phase-shift between 
pairs £ A-B’ and £ C-D\ 


5.2. Relation to BDW in cuprates 

Our aim is to investigate the nature of density-wave instabilities for the cuprate-like 
Fermi surface. A number of recent works 41,47,97 99 have analyzed these instabilities 
for the full Fermi surface (Fig. 16a), starting with the work in Ref. 41, which analyzed 
the form-factors and wavevectors of the (sub-)leading instabilities within a Hartree- 
Fock computation for H = Ht + Hj. The lowest eigenvalue, Aq, of the particle-hole 
T-matrix determines the structure of the charge-order and is shown in Fig. 16 as a 
function of Q. 

Remarkably, the leading density-wave instability occurs for a wavevector Q b = 
(Q o, ±Qo) 5 as shown in Fig. 16a. The form-factor for this BDW is purely d-wave, as 
required by symmetry. Moreover, Q b tracks the diagonal hot-spot wavevector over 
















January 5, 2015 1:25 


WSPC Proceedings - 9.75in x 6.5in 


DCSS'arxiv page 27 


27 


a range of doping. Thus even for the lattice Hamiltonian, and far away from any 
putative AFM QCP, the leading instability in the particle-hole channel corresponds 
to the instability that we had obtained earlier in the hot-spot theory. The one 
important difference however is that this instability is not degenerate with the d- 
wave superconducting instability, as expected in the absence of the exact SU(2) 
symmetry. The eigenvalues of the T-matrix in the particle-particle channel were 
also analyzed 100 , allowing for the possibility of a finite center-of-mass momentum, Q 
(i.e. a pair-density wave, PDW, state 101,102 ). However, the leading weak-coupling 
instability was always to a d-wave SC with Q = 0. It has been argued recently that 
there could be a strong-coupling instability to a PDW state, without reference to 
any microscopic model, arising from an effect analogous to ‘Ampere’s law’ 103 . 

Coming back to our discussion of Aq, one can also notice the ‘valleys’ of local 
minima that extend from Q b to Q x a = (Qo, 0) and Q y a = (0, Qq). The form-factor for 
these BDWs predominantly have a d-wave component with a small 5, s'— component. 
Notice that the Q x,y are oriented along the Cu-0 directions and are therefore quite 
similar, at least qualitatively, to the states observed in the experiments. However, 
while here the wavevectors connect the hot-spots, it is unclear to what extent this 
is the case experimentally across different families of the cuprates. 





- 0.080 


- 0.274 


Fig. 16. Left: A cuprate-like Fermi surface with t\ = 1.0, t<i — —0.32, £3 = 0.128, \i — —1.11856. 
Black circles represent the hot-spots. The arrows represent wavevectors Q ® = (Qo,0) and Qb = 
(Q 0 , Qo) corresponding to BDW-a and BDW-b. Right: A q as a function of Q, showing the global 
minimum at Q b and two ‘valleys’ of local minima extending from Q b to Q®, Qa • Another local 
minimum appears close to Q = ( 7 r, 7 r), which also breaks time-reversal symmetry and represents 
the ‘staggered-flux’ state 5 . Adapted from Ref. 41 

The above computation was somewhat rudimentary—it only analyzed the 
quadratic instabilities of the t-J model (without projection). A subsequent study" 
carried out a variational Monte-Carlo computation for the leading instabilities of 
the t-J model, with an infinite on-site Hubbard U (implemented using Gutzwiller 
projection). Their results are shown in Fig. 17. Interestingly, there were regimes in 
parameter space where for at least some Fermi surface geometries, the BDW state 
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with Q a was found to be the leading instability. However, this required somewhat 
large nearest neighbor Coulomb repulsion (V J); furthermore, given the large 
coupling and the numerical projection, we cannot be certain that this instability is 
arising from a Fermi liquid with a large Fermi surface. 


+ 0.00 + 0.00 + 0.00 



Fig. 17. Top row: Gain in variational energy for ordering at wave vector Q; few representative 
parameter sets have been shown, corresponding to axial BDW, staggered flux and diagonal BDW 
respectively. The red dots denote the hot-spot wavevectors. Bottom row: Approximate diagrams 
of the leading instability in the particle-hole channel, as a function of J and V. Adapted from 
Ref. 99. 

The above analysis leads to two key insights, which is responsible for the diag¬ 
onal BDW with a d-form factor being the leading density-wave instability: ( i ) the 
remnant of the SU(2) pseudospin symmetry, arising from the exchange interactions, 
after it has been broken explicitly is still strong enough, and, (ii) a large joint den¬ 
sity of states is available for hot-spots along the diagonal wave vector, which favors 
the diagonal BDW. One possible way around this is to preferentially suppress the 
density of states in the antinodal regions. 

A natural ingredient that suppresses density of states at the antinodes is su¬ 
perconductivity. However, since the BDW correlations onset above T c , one can 
only allude to superconducting fluctuations. Indeed, pairing fluctuations in a metal 
suppress the diagonal BDW much more strongly compared to the axial BDW at 
low enough temperatures 100,104 . On the other hand, experiments have detected the 
presence of these pairing fluctuations 16 to a relatively narrow window above T c , 
while the onset of BDW occurs significantly above T c . Moreover, this approach still 
does not rule out why the diagonal BDW has never been detected in the cuprates 
at high enough temperatures. 
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Therefore, one needs a different mechanism to suppress the antinodal density of 
states. One other possibility is to use long-range commensurate AFM to gap out the 
antinodes completely 45,46 . The instabilities of a reconstructed Fermi surface with 
hole-pockets (Fig. 10) does give rise to a density-wave with the correct wavevec- 
tors, but with a form-factor that is predominantly s /45 . However, this can not be 
an explanation for the charge-ordering phenomenon in the non-La-based cuprates 
where the AFM correlation length is less than one lattice constant. 

It is then perhaps appropriate to investigate the instabilities of a ‘normal’ state 
that already has a well-developed pseudogap, without any symmetry-breaking, i.e. 
one should investigate the instabilities of a state that has Fermi-‘arcs’, and not a 
large Fermi surface. This naturally leads us to investigate the instabilities of the 
FL* introduced in Section 4. 


5.3. Ordering instabilities of FL* 

Computing the ordering instabilities of the FL* metal is a far more complicated 
task than that for the Fermi liquid described in Section 5.1. Apart from instabilities 
associated with particle-particle and particle-hole excitations near the Fermi surface, 
we also have to consider the instabilities associated with the gauge sector. For the 
U(1)-FL* case, the latter contain monopole tunneling events which necessarily lead 
to confinement and broken symmetries at low temperatures. There has been partial 
progress on these difficult issues in the literature, but we shall not review them here. 
Instead we shall focus only on the charge ordering instabilities associated with the 
‘small’ Fermi surface, in close analogy to the analyses of Section 5.1. 

The surprising outcome is that even a somewhat rudimentary RPA level analysis 
of the density-wave instabilities of the FL* immediately leads to a state with d-form 
factor and wavevectors along the axial directions, as displayed in Fig. 18. Moreover, 
the wavevector is related to a geometric property of the underlying Fermi surface 
— it is determined by the distance between the tips of the pockets. Hence, as the 
hole-doping increases, the magnitude of the wavevector decreases, consistent with 
the experimental results on the non-La-based cuprates 105 . 

The above discussion seems to suggest that the d-form factor BDW with the 
axial wavevectors, observed experimentally at low temperatures, arises most natu¬ 
rally out of a parent high-temperature FL* phase. However, the plain vanilla FL* 
state has hole pockets with nonzero (but small) spectral weight on the backside. 
While the photoemission results of Ref. 10 have been argued to be approximately 
consistent with this, clear signatures of the backside continue to be elusive. A 
promising future direction for experiments should be to employ probes which do 
not involve adding or removing an electron from the sample (which are sensitive to 
the quasiparticle residue), but instead measure the 2k p response corresponding to 
the ‘small’ Fermi surfaces directly using Friedel oscillations, the Kohn anomaly, or 
ultrasonic attenuation. 

In Section 4, both prescriptions of arriving at the FL* relied on starting from the 










January 5, 2015 1:25 


WSPC Proceedings - 9.75in x 6.5in 


DCSS'arxiv page 30 


30 



Fig. 18. Left: The spectral function for a FL* state with the same FL parameters as in Fig. 16. 
The other parameters are to = —0.5ti, t\ = 0.6ti, A = 0.75ti (Eq.27). The black arrows denote 
the wave vectors (Q a ) of the leading density-wave instability. Right: A q as a function of Q, 
showing the global minimum at Q a . Adapted from Ref. 47. 

antiferromagnetic insulator and doping it with charge-carriers. A natural question 
that now needs to be addressed is how to connect the metallic FL* phase with a 
‘small’ Fermi surface to the ‘large’ Fermi surface metal. As discussed recently 82,84 , 
it is possible to describe this in terms of a Higgs’ transition, that does not involve 
any broken symmetries and yet reconstructs the Fermi surface. However the tran¬ 
sition proceeds via an intermediate non Fermi-liquid phase which was dubbed an 
SU(2) algebraic charge-liquid (ACL). Our view of the proposed phase diagram is 
reproduced in Fig. 19 for the sake of completeness. We refer the interested reader 
to Ref. 82 for a more complete discussion of the Higgs’ transition and a comparison 
to some of the associated phenomenology in the cuprates. 

6. Outlook 

The main purpose of this review was to address a number of pertinent questions 
regarding the nature of the newly discovered charge-order in the (non-La-based) 
underdoped cuprates and to explore its connection with the high temperature pseu¬ 
dogap phase. Let us revisit the questions that we had raised in the introduction 
one by one, in order to summarize our current understanding and provide a brief 
outlook: 

• The d-form factor BDW likely arises as a low-temperature instability of 
a fractionalized Fermi-liquid (FL*) with ‘small’ Fermi surfaces of hole-like 
quasiparticles 47 . While the FL* resembles a usual Fermi-liquid in terms of 
its transport and thermodynamic properties, it violates Luttinger’s theorem 
on the area enclosed by the Fermi surface, and this is closely linked to the 
presence of emergent excitations belonging to a gauge sector 80 . 

• The FL* has pocket Fermi surfaces with low intensity on the ‘back sides’ be- 
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T 



QCP 


Fig. 19. Our proposed global phase diagram for the cuprates, adapted from Ref. 82. The present 
review has focused on the regime to the left of the Higgs quantum critical point (QCP), with the 
U(l) FL* and its instability to the d-BDW. The Higgs QCP describes a phase transition from 
the U(l) ACL to the SU(2) ACL which is accompanied by a change from ‘small’ to ‘large’ Fermi 
surfaces. 

cause of small (but non-zero) quasiparticle residue 72,87 . This is presumably 
the reason only ‘Fermi arcs’ have been detected by photoemission so far. 
Fermi surface probes which don’t involve adding or removing an electron 
could provide detection of the full pockets. 

• The modulation wave vector of the BDW is along the axial directions and 
is determined by the distance between the ‘tips’ of the pockets 106,107 . The 
d-form factor arises as a result of the remnant SU(2) particle-hole symmetry 
associated with the exchange interactions 42 , that maps d-wave supercon¬ 
ductivity to d-form factor BDW. 

• The non-mean field nature of the onset of BDW and SC correlations is 
governed by the angular fluctuations of a six-component order-parameter. 
A phenomenological non-linear cr-model for the multicomponent order- 
parameter provides a quantitative description for the onset of charge-order 
intensity as well as the diamagnetic response in the pseudogap regime of 
the underdoped cuprates 57,60 . 

• A large Fermi surface in the presence of long-range (and nearly bidirec¬ 
tional) charge order leads to a nodal electron pocket 62 flanked by recently 
predicted small hole pockets 61 . These are likely responsible for the observed 
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quantum oscillations 63 . At low temperatures and high fields, we expect the 
U(l) FL* to undergo a crossover into a confined phase with charge order; 
it is possible that the nodal Fermi surface region is insensitive to confine¬ 
ment and so has a pocket structure similar to that obtained from the large 
Fermi surface computation in Ref. 61. It is clear that this crossover to 
confinement is an important topic for future research. 

• A candidate theory 82 for the optimal-doping QCP is associated with a 
Higgs’ transition, without any broken symmetries, from a SU(2) ACL to a 
U(l) ACL. 
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